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Abstract 

We calculate the entanglement between two spins in the ferromagnetic Heisen- 
berg chain at low temperatures, and show that when only the ground state and 
the one particle states are populated, the entanglement profile is a gaussian with 
a characteristic length depending on the temperature and the coupling between 
spins. The magnetic field only affects the amplitude of the profile and not its 
characteristic length. 
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1 Introduction 



Thermal entanglement in many body quantum systems has recently attracted a lot of 
attention ^[21IHllUli3iniEllBllHl- This term is usually used to specify the amount of 
entanglement which exists between two spins of a multi spin system when the whole 
system is in a state of thermal equilibrium. There are at least two very good reasons 
for the interest in this problem. First, it is well known that many body quantum 
systems may be more correlated than the corresponding classical systems, the excess 
correlation stemming from a very basic quantum property, namely quantum entan- 
glement |1U1 1111 112) . Thanks to the progress in the past few years, we have now 
good measures of entanglement (at least in certain limited cases ^Hl) an d it is most 
natural to quantify this excess correlation or entanglement in many body systems. 
Once a property like entanglement becomes quantifiable, we can use it as a kind of 
order parameter and study its behavior when the external parameters of the system 
change. This will certainly shed light on the nature of quantum phase transition 
(qualitative change in the nature of the ground state) ^1] . The second reason is that 
one dimensional arrays of spins, are a natural candidate for storing quantum infor- 
mation |15l I16j . Entanglement is a valuable resource for implementation of any kind 
of quantum algorithm and quantum information protocol |17( I18j and one needs to 
find the amount of entanglement which exists between spins or qubits and the ways 
to control it, when the array is in thermal equilibrium. 

A prototype of a many body system is the isotropic Heisenberg spin chain whose 
hamiltonian is as follows: 

N-i N-l 
H = -Jj2^r + ^ E °i ■ (!) 

1=0 1=0 

Here af, where (a = x,y or z) is a spin (Pauli) operator acting only on site I and 
periodic boundary conditions are assumed. 

In this context and to our knowledge, up to now the following studies have been 
made. Nielsen |3] has studied the entanglement of a two-spin system at finite tem- 
perature interacting via the Heisenberg interaction. Wang 0IEj and Rigolin 19 have 
also studied the effect of anisotropy on the thermal entanglement in a two spin sys- 
tem. Arnesen, Bose and Vedral have numerically studied the variation of two spin 
entanglement with temperature and magnetic field in a few-spin isotropic Heisenberg 
chain and Gunlycke et al. |2j have studied the same problem in an Ising model in a 
transverse field where they have observed a kind of quantum phase transition at zero 
temperature when the entanglement suddenly shows up even with an infinitesimal 
transverse magnetic field. O'Connor and Wootters (2131) have calculated the maxi- 
mum possible entanglement between nearest neighbors in a translationaly invariant 
state and have shown that the ground state of the antiferromagnetic Heisenberg chain 
satisfies this maximum under certain conditions. Osterloh et al [7j have shown that 
quantum phase transitions of a class of spin systems can be characterized by the 
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change in the entanglement between the next and next-nearest neighbors spins. In 
particular they have shown that near the point of quantum phase transition the 
nearest and next-nearest neighbor entanglement exhibit logarithmic divergence and 
universal behavior. A similar study has been done in j^j. 

Of particular interest to us here is the work of Wang and Zanardi jlj who have 
shown that in the isotropic Heisenberg model in the absence of magnetic field, the 
nearest-neighbor entanglement can be related directly to the free energy of the model. 
Thus by knowing only the eigenvalues and not the eigenstates, one can calculate the 
entanglement between nearest neighbor spins. In particular they have shown that 
in an Heisenberg ferromagnet there is no entanglement between nearest neighbors if 
B = 0. 

This is a very rare and fortunate situation where one can draw general conclusions 
about a problem whose solution generally needs a knowledge of the whole spectrum, 
i.e. the eigenvalues and eigenvectors. As they have correctly pointed out, their argu- 
ment is based on translational invariance of the lattice and the su(2) symmetry of the 
model. Once the su(2) symmetry is broken, e.g. by applying a magnetic field, their 
argument will not be valid anymore and one may find entanglement in the Heisenberg 
ferromagnetic chain. 

In this paper we show that if we apply a magnetic field to a ferromagnetic Heisenberg 
chain, then pairwise entanglement will develop between spins at arbitrary given sites. 
An exact solution of this problem is extremely difficult if not impossible, since as we 
will see it requires the determination of the spin-spin correlation functions over all 
the energy eigenstates . However at low temperatures when only the ground state 
and the first excited states are populated one can calculate the entanglement ana- 
lytically. Moreover one can now calculate the entanglement profile (entanglement 
between arbitrary sites) and see its dependence on various control parameters, like 
the spin coupling, the magnetic field and the temperature. 

We will find that for finite but arbitrary number of spins, the entanglement profile is 
a guassian with a characteristic length depending only on the temperature and the 
coupling between spins. On the other hand the magnetic field and the total number 
of spins only affect the amplitude of the profile. 

Remark Contrary to the case where there is no magnetic field, an exact proof 
on the presence or absence of entanglement in the Heisenberg ferromagnetic chain is 
still missing. Exact solutions for up to 4 and numerical evidence for up to 10 spins 

indicate that there is no entanglement in ferromagnetic chains. If this is indeed 
the case then the entanglement calculated in this paper should be looked upon not 
as a result of approximation but as a result of truncation of the spectrum. That 
is if we somehow prevent the higher excited states from mixing with the low lying 
states of the Hiesenberg ferromagnet, then we can produce entanglement between 
different remote sites and suitably manage this entanglement. This problem may be 
of practical relevance in solid state implementations of quantum computers which 
should be kept at low temperature to reduce the effect of noise. In such cases we need 
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to know not only the entanglement between nearest neighbor spins but also between 
distant spins. The variation of entanglement with distance has not been studied in 
the works mentioned above, due to the complications in the ground and excited states 
of the anti-ferromagnetic Heisenberg chain or the Ising model in transverse field. The 
low lying states of the Heisenberg ferromagnet is a simple situation where we can 
study this variation. 

The structure of this paper is as follows: In section |2] we discuss the general form of 
the two body density matrix in the Heisenberg chain. In section |3] we review the basic 
properties of the low lying states of the Heisenberg ferromagnet and in section |3] we 
calculate the explicit form of the two particle density matrix and its concurrence which 
characterizes its degree of entanglement. We conclude the paper with a discussion in 
section |5J 

2 Some general properties of the two particle density 
matrix 

The Hamiltonian of the Heisenberg spin chain (^Q) has the following symmetries 

[H,S z ] = [H,tt] = [H,T]=0, (2) 

where 

1 N - 1 

S z ■= ? E °f. (3) 

is the third component of the total spin, and T is the translation operator 

T\s , si,--- s N -i) = \si, s 2 , ■ ■ ■ sjv-i, s ). (4) 

Here |sj), s, = ±1 is the state of the i-th spin expressed in terms of the eigenvectors 
of a z '', i.e. of |sj) = Sj|sj) and tt is the operator which reflects the spin chain around a 
suitable site: i.e. 

tt\s , si,--- S N -l) = \s N -i,s N - 2 , ■ ■ ■ s ). (5) 

This symmetry means that enumerating the sites of the lattice in the clockwise or anti- 
clockwise directions is immaterial for the derivation of the properties of the lattice. 
The equilibrium state of such a system at temperature T is given by a density matrix 
p = 6 z , where = jr^f, &b is the Boltzman constant and Z = tr(e~^ H ) is the 
partition function. 

We are interested in the reduced density matrix and entanglement of two spins at 
sites m and n. This is given by 

Pm,n = ^tr^e-P", (6) 

where fn^ri means that the trace is taken over all factors in the tensor product space 
except those at sites m and n. 

From © we find various elements of the matrix Pm^n &s 
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(8) 



1 

(ij\p m ,n\k,l) = -tr (e~P H \E u {m)E lj {ri))) , (7) 

where i, j, fc, Z = 0,1 and Eki(m) is the operator which acts like E^i = \k){i\ on site m 
and like identity on all the other sites. (Note that we use two equivalent notations for 
the states, namely |+) = |0) and |— ) = |1), the ± notation is usual for the description 
of spin states and the |0) , 1 1) notation is usual in description of qubit states in quantum 
computing.) Expressing the matrices E^ in terms of Pauli matrices and using the 
symmetry [H, S z ] = 0, one can rewrite p m . n in terms of correlation functions of spin 
chains as: 

/ u+ \ 

w z 
z w 

V u ~ 1 

where due to translational and reflection invariance, the parameters u^, w and z 
depend on \m — n\ and 

u- := ^tr((l-^)(l-<)e-^) 
w := ifr((l-<)(l + ^"' W ) 

z ■= ^tr(a m a+e-P H ). (9) 

For notational convenience we do not show the explicit dependence on \m — n\ in 
the parameters of p m , n - The equality of the elements /5oi,io and Piofli, (the reality of 
z) is a consequence of the tt symmetry and translational symmetry. This is pointed 
out by Wang and Zanardi in 4 , although they explicitly derive this property by the 
symmetry {(a x )® N , H] = which exists only in the absence of magnetic field. To see 
this we temporarily insert the site dependence in the parameters and note that 

z(m -n) = tr(pcr-a+) = tr^pTrcr^a^) 

= tr(pira m a+7r) = tr(po N _ m o^_ n ) 

= tr(pa+_ n a„_J = z*(N-n-(N-m)) 

= z*(m-n). (10) 

The same argument shows why poi,oi = Pio,io- 

The eigenvalues of this density matrix (which is positive) are , and w ± z, which 
implies that the u^ 1 and w are positive, while the sign of z is not determined. 
The entanglement of the two spins at sites m and n is determined by the concurrence 
C(|m — n|) defined as 

C(\m — n\) = max{0, X\ — A2 — A3 — A4}, (11) 

where Ai, A2, A3, and A4 are the positive square roots of the eigenvalues of the matrix 
Pm,nPm,n in decreasing order. The matrix p is defined as 

Pm,n = {<r y (8) a y )p*(a y ® a y ), (12) 
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where * denotes complex conjugation in the computational basis. It turns out that 
[2~U] C(\m — n|) is given by 

C(\m - n\) = 2 max (0, |z| - Vu+u"). (13) 

3 The low lying states of the Heisenberg ferromagnet 

The Hamiltonian Q can be rewritten as 

N-l N-l 
H = JN-2JY / Pl.l+1 + E °u (14) 

1=0 1=0 

where Pij+i is the permutation operator acting on sites I and I + 1, P\s, s') = \s', s). 
For the ferromagnetic chain (J > 0), the ground state of H is the disentangled state 

|eo) = |-, ',-,-,-), (15) 
where all the spins are down, with energy eo = — (J + fiB)N. 

The first excited states are linear combination of one particle states 

iv>> = E c ^>> (i 6 ) 

k 

where 

,-,-->, (17) 
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in which only the spin in the k— th site is up. 
From (|14|) we find 

£T|ife) = -2J|jfe - 1) - 2J\k + 1) + (J(4 — JV) + /iS(2 - N)) \k). (18) 

leading to the recursion relations for the coefficients 

- 2J( Cfe _! + c k+1 ) + (J(4 - N) + M £(2 - AT)) Cfe = Ac fc (19) 

whose solution together with the periodic boundary condition gives us the final form 
of the first excited states labeled by \ip s ) with s = • • • N — 1: 

N-l 

i^) = - 7 ^E e ~i A; )' ( 2 °) 



o 

with energies 

\ s = eo + 2fiB + 8Jsin 2 ^. (21) 

The lowest energy in this band belongs to the flat wave (s = 0) which is higher than 
the ground state energy by 2fiB. In the classical picture this is the energy difference 
due to the flipping of one spin in the magnetic field and there is no contribution from 
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the interaction of spins. One can guess that the lowest energy in the two particle 
sector also corresponds to a flat wave, namely to a state 



I*) 




0<i<j<N-l 



(22) 



where \i,j) is a state in which only the spins at sites i and j are up. It is easy to see 
that this is indeed an eigenstate of the Hamiltonian with an energy eo + AfxB. 
Thus if ^ » 1, one may assume that only the ground state and the first excited 
states will be populated. We will assume that the magnetic field is strong enough 
and the temperature is low enough so that this condition is fulfilled. 



4 Thermal entanglement of two spins 

We will calculate the density matrix © at low temperatures by taking into account 
only the ground state and the one-particle eigenstates. At low temperature we ap- 
proximate various quantities in Q as follows: 

N-l 

tr(Ae^ H ) « (eo\A\eo)e-fa + £ (V^l^e"^. (23) 

s=0 

The operator A takes one of the following forms, (l±0^j(l±erjj), — a*), 

or °mO 

We need the relevant matrix elements over the one particle states. It is easily verified 
that 



~(I + OI*) = 5 n , k \k), 

l -(I-a z n )\k) = (l-6 n , k )\k), 

a-\k) = 6 n;k \k), (24) 

from which we obtain the following matrix elements, where we have used the fact that 
m 7^ n: 

1(^1(1 + 0(1 + 01^) = 0, 

1 1 N — l 2-Ki{l — k)s \ 

-(^1(1+0(1-01^) = -T7 E e~^~6m,k8k,l = T7, 

k,l=Q 

1(^1(1-0(1-01^) = 1-1 

, I -\-2iri(m — n)s 

(tp s \a-aZ\^ s ) = j^e » . (25) 
From Q and (|25|) we find values of different matrix elements of p mn as follows: 
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u + = 

1 El-o 1 



w = 



N 



1 1 JV-i 

z = ±±y e — SH~^«. (26) 
The partition function at this level of approximation is given by 

N-l N-l 
Z « e -^o + £ e -/3A s = e -/3,o + e -/?(6o+2^) £ e -8/Wsin» (2?) 

fc=0 s=0 

For large iV we can approximate the sums in the above formulas by integrals 
and at low temperatures (large (3) we can evaluate the integrals by saddle point 
approximation. Thus 



N-l N-l 

k=0 s=0 

„ e -/3(,o+2MB) / f e -8/3Jsin^ ds 

« e-^+^ivJA- . (28) 

The partition function will be 

(l+^tf-^). (29) 

We also find 

AT— 1 JV 

e _ 2 li (n- Z n) £ _ /3A5 ^ e _ /3(eo+2/ , B) CI & _^J gin 2 ™ _ -2,i(n-m). ^ 



s=0 ,/ ~ 

e 



V 8vr/3J v ; 



Putting all this together we find 

/ o 



1 

Pm,n 



1 e 8 /3 J 

(re — m)^ 



e 8 /" 1 
V iV - 2 + ^/Wj5Je 2 ^ B J 



(31) 
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Now that we have the two particle density matrices, we can calculate the entan- 
glement of two spins at an arbitrary distance. From Q13JI we find 



C(\m-n\) = 2\z 



2 



e 




(32) 



N + V8¥ple 2 ^ B 



which can be rewritten as 



C(\m-n\) =C e-i 



1 (m — n) 2 
2—P 



(33) 



where 



C : 



2 



and 




(34) 



denote respectively a scale of the value of entanglement and a kind of entanglement 
length which determines a scale of length over which the distant spins are still entan- 
gled. It is seen that the magnetic field has no effect on the entanglement length and 
only controls the entanglement amplitude. The higher the magnetic field the lower 
the value of Co . This is understandable since a high value of the magnetic field orders 
all the spin more in the direction of the z axis in favor of the disentangled ground 
state. On the other hand it is the coupling between adjacent spins that entangles 
distant spins, hence the entanglement length is only affected by J and T. Raising 
the temperature has a destructive effect on the entanglement length since thermal 
fluctuations tend to move the system away from the quantum regime. Finally we 
note that in the thermodynamic limit when N — ► oo, the concurrence Cq tends to 
zero, so that even the nearest neighbor sites will not be entangled anymore. Figure 
^ shows the entanglement profile for a spin chain of N = 20 spins at two different 
temperatures. 

5 Discussion 

As Wang and Zanardi [I] have shown, in the Heisenberg ferromagnet in the absence 
of a magnetic field there can be no thermal entanglement between nearest neighbors 
spins. Their argument is based on the su(2) symmetry of the Hamiltonian. In the 
presence of a magnetic field this symmetry breaks down to a u(l) symmetry and one 
can only derive entanglement by explicitly calculating the relevant thermal correlation 
functions of the model. We have shown that at low temperatures or high magnetic 
fields where to a good approximation only the ground state and the one particle states 
are populated, entanglement can exist between different sites. We have shown that 
the entanglement profile is a gaussian with an amplitude depending on the magnetic 
field and temperature and a characteristic length depending only on temperature and 
the coupling between spins. 

Another important point is that in the previous works the pairwise entanglement has 
been studied in more complicated situations and for making connections to quantum 
phase transitions, but only for nearest neighbors. By considering the forromagnetic 
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Figure 1: The concurrence C as a function of x = \m — n\, the distance between 
spins, for a spin chain of iV = 20 spins at two different temperatures. Solid line 

(w = 3 > w = °- 6 )' dashed line (w= 4 >w = °- 8 ) 
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case and by restricting to the low lying states, we have been able to see how en- 
tanglement varies with distance. Knowing this variation may be important for the 
implementation of quantum protocols on arrays of qubits in any candidate for imple- 
mentation of quantum computers. An analogous study even for the ground state of 
the anti- ferromagnetic spin chain is extremely difficult, since it requires the calcula- 
tion of correlation functions for arbitrary distances and up to now the exact ground 
state correlation functions of the anti-ferromagnetic model have been calculated for 
distances of only up to 4 sites [22] . 
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